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ABSTRACT 

We demonstrate the unique capabilities of Herschel to study very young luminous extragalactic young stellar objects (YSOs) by analyzing a 
central strip of the Large Magellanic Cloud obtained through the HERITAGE Science Demonstration Program. We combine PACS 100 and 160, 
and SPIRE 250, 350, and 500 yt/m photometry with 2MASS (1 .25-2. 17 jdm) and Spitzer IRAC and MIPS (3.6-70 jdm) to construct complete spectral 
energy distributions (SEDs) of compact sources. From these, we identify 207 candidate embedded YSOs in the observed region, ~40% never- 
before identified. We discuss their position in far-infrared color-magnitude space, comparing with previously studied, spectroscopically confirmed 
YSOs and maser emission. All have red colors indicating massive cool envelopes and great youth. We analyze four example YSOs, determining 
their physical properties by fitting their SEDs with radiative transfer models. Fitting full SEDs including the Herschel data requires us to increase 
the size and mass of envelopes included in the models. This implies higher accretion rates (>10~'^Moyr~^), in agreement with previous outflow 
studies of high-mass protostars. Our results show that Herschel provides reliable longwave SEDs of large samples of high-mass YSOs; discovers 
the youngest YSOs whose SEDs peak in Herschel bands; and constrains the physical properties and evolutionary stages of YSOs more precisely 
than was previously possible. 
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1. Introduction 

The proximity of the Magellanic Clouds offers a unique op- 
portunity to analyze the complete inventory of luminous YSOs 
over an entire galaxy. With known YSO distances, luminosi- 
ties, masses, and mass accretion rates can all be well-defined. 
Comparison of the properties of YSOs in the Magellanic Clouds 
and in the Milky Way can reveal diflferences in star formation 
physics due to metallicity and environment. 

Using the Spitzer SAGE ("Surveying the Agents of Galaxy 
Evolution") Survey of the Large Magellanic Cloud (LMC; 
iMeixner et alll2006k . Whitney et al. (2008; W08) and Gruendl & 
Chu (2009; GC09) discovered -1800 massive YSO candidates 
in the LMC (a 90-fold increase over previous work). Spitzer 
studies selected sources using colors and SEDs at wavelengths 
<24 yum (where Spitzer can resolve individual YSOs), requiring 
a detection at 4.5 yum or shorter in most cases. These surveys 
thus missed the youngest, most embedded YSOs that can only 
be detected at longer wavelengths. 

The Herschel Space Observatory ( Pilbratt et al.l '2010) has 
the spatial resolution required to study individual sources at 
X >50 yum (from - 1.3 pc at 70 // m to -8.7 pc at 500 //m for 
a distance of 50 kpc, Schaefer 2008). The least-evolved massive 
protostars are characterized by cold dust temperatures probed at 
far-infrared (far-IR) wavelengths, and are expected to be -10^ 
times brighter at 100 yum than at 5 yum (IWhitnev et al.l 120041: 
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iMolinari et al.l2008l ). making Herschel extremely eflTective at de- 
tecting those youngest YSOs. 

With Herschel, we not only discover new objects but 
also better characterize 5/7/^z^r-identified YSO candidates. The 
Herschel data constrain the models of these sources and im- 
prove estimates of such physical parameters as total luminos- 
ity, stellar mass, and total dust mass. We demonstrate these ca- 
pabilities by studying a strip across the LMC observed as part 
of the Science Demonstration Program (SDP) - the first part of 
the Herschel Key Program "HERsc hel Inventory of the Agents 
of Galaxy Evolution" (HERITAGE: iMeixner et al. 1'2010') in the 
Magellanic Clo uds. The strip was ma pped in the PACS 100 and 
160 yum bands (Poglitsch et al. 2010) and SPIRE 250, 350, and 
500 yum bands (.Griflan et al...2010.) . 

2. Source selection and photometry 

An initial list of 640 sources was compiled by hand, choos- 
ing apparent point sources in Herschel images. Astrometry was 
refined using the SAGE-LMC MIPS 24 yum image. We per- 
formed aperture photometry on 2MASS, Spitzer, and Herschel 
images using apertures scaled to the instrumental resolution: 
A" radius for 2MASS, y for IRAC, 6'' and 12'' for MIPS 
24 and 70 yum, 8'' and 12'' for PACS 100 and 160, and 13", 
17", and 23" for SPIRE 250, 350, and 500, respectively. The 
aperture size relative to PSF was chosen after examination of 
the multi-wavelength images, to most consistently measure a 
single YSO as distinct from its environment. Aperture correc- 
tions are 1.4, 1.5, 1.7, and 1.7 for IRAC, 1.8 for MIPS, 1.4 
for PACS, and 1.3, 1.25, and 1.2 for SPIRE, which not only 
account for the portion of the PSF extending outside the aper- 
ture, but also flux contamination from the low-level PSF wing 
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in the background annulus. Background emission was calcu- 
lated as the sigma-clipped mean of an annulus spanning 1.75 
to 2 times the source radius. Flux uncertainties were increased 
where large gradients existed across the background annulus and 
in regions of crowding and confusion. Photometry was verified 
quantitatively through comparison to SAGE PSF-fit photome- 
try and by manually measuring -20 sources (including those 
whose SEDs are highlighted below). No significant systematic 
off'sets were found, and random off'sets are consistent with the 
quoted uncertainties. We find that in these early data we can re- 
liably extract point sources as faint as (200,300,150,70,40) mJy 
at (100,160,250,350,500) yum. Greater integration time and im- 
proved artifact mitigation will result in fainter values for the 
complete HERITAGE survey. T he uncertainty in the absolute 
flux calibra tion is 20% for PACS (iPoglitsch et al.l2010l) and 15% 
for SPIRE (ISwinvardetal.1120101) . 

We selected a subset of reliable Herschel sources, most likely 
YSO candidates, from the more complete source list using the 
following criteria. We carefully examined the environment of all 
candidates at all wavelengths, simultaneously with their SEDs, 
aperture photometry, and existing W08 and GC09 catalog pho- 
tometry. Sources were removed that could not be unambigu- 
ously identified in images over a wide wavelength range due to 
multiplicity or complex diff'use emission. We also required re- 
liable photometry spanning at least 5.8-100 yum. This examina- 
tion is subjective but consistent with the goals of assessing how 
Herschel changes our understanding of star formation. Analysis 
with quantitative completeness limits will be performed when 
the higher quality HERITAGE survey data are available. 

We removed known non-YSOs from our list. At Herschel 
wavelengths, background galaxies are the main contaminant as 
their far-IR emission arises fro m star-forming regions. We ex - 
cluded 2 AGN candidates from lKozlowski & KochaneS (l2009h . 
as well as 7 probable galaxies based on spatial morphology in 
high resolution IRAC images. Evolved stellar envelopes lack 
cold dust, peak in the mid-I R, and are few amo ng Herschel 
sources. We removed one from IB oyer et al.l (l2010h . leaving 207 
sources in our very conservative list. 



3. Far-IR photometric properties of LMC YSOs 

Figure [T] shows a color-magnitude diagram of our 207 YSO 
candidates, combining Spitzer (70 yum) with Herschel (160 
and 250 yum) bands. Other color combinations show simi- 
lar source distributions and lead to similar conclusions. YSO 
candidates selected based on photometric data alone require 
spectroscopic confirmation. The presence of ice absorption is 
quite definitive; silicate absorption and aromatic emission are 
at least strongly suggestive (in the parsec- sized beam). To 
date, -300 sources in the LMC have been spectrosco pically 
confir med or supported jShimonishi et al. 2008; Oliveira et al.l 
I2Q09I: ISeale et al.lf2009l: R^n Loon et al.f 2010). Of these, 126 lie 
in the HERITAGE SDP area; our list includes 58 (the remainder 
are too faint or or too confused). We label these as "confirmed 
YSOs." Photometric YSO candidates may also be "confirmed" 
through maser identification. In the HERITAGE SDP area, 
maser emission has been detected in the N 113 (H2O and OH; 
Whiteo^¥G^dneilfT98F 
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„ iLazendic et al.ll2002t lOliveira et"aD 

2006; Brooks & Whiteo ^ll997l see Fig H and N 105 (H2O 
OH and methanol; iSinclair et al.l Il992t ILazendic et al. I l2002t 
lOliveira et al.l l2006l) star formation regions. Each of these re- 
gions has two clumps of masers, and all four sites are associ- 
ated with mid-IR sources identified as YSO candidates. Three of 
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Fig. 1. Color-magnitude diagram showing Fpacs miumfpMiPSiOfim 
vs. FspiRE250i^m Uy]- The distribution of YSOs is compared to 
predictions from YSO models shown in grey (R06); those to the 
right of the dashed line are Stage O/I models. We also compare 
to optically thin greybody emission with M^^^ = (0.5-150) x 
10^ M© at 20-40 K, showing that these sources have very cool 
and massive circumstellar envelopes. The Herschel sources dis- 
cussed in Sect. [5] are indicated by red diamonds. A typical error 
bar is shown in the lower left corner. 



these sources have been spectroscopically observed (Scale et al. 
2009). 

We fit the SED of each source with a large grid of 
dust radiative transfer mo dels of individual massive YSOs 
(iRobitaille et al.ll2006L 12007 . hereafter R06, R07), constraining 
the circumstellar dust distribution and other physical parameters 
(Sect.O. Comparison between observed far-IR fluxes and those 
predicted by the R06 YSO model grid confirms the youthful na- 
ture of these sources. All Herschel sources, and all models to the 
right (redward) of the dashed line in Fig. [T] are consistent with 
Stage O/I , i.e. the circumstellar mass exceeds that of the central 
source. In fact, many measured fluxes are redder and brighter 
than the R06 grid. We model example sources in detail below, 
and in general this color difl'erence is resolved by adding a larger 
outer envelope of cool dust. 

We also estimate the mass of cool circumstellar dust using 
simple greybody emission from optically thin dust at a single 
temperature. A temperature - mass grid is overplotted on Fig.[T] 
We use the same opacity curve as the YSO radiative transfer 
models, i.e. the opacity power-law index J3 is interstellar (Fj o c 
A'^ withyS ^ 2) for wavelengths > 20 yum (Whitney et al.l2003bh . 



4, Herschel view of YSOs and their environments 

We highlight four YSOs, showing the difl'erent ways that 
Herschel can impact star formation studies (SEDs are modeled 
in Sect. O. Two sources are in the N 113 Hii region, one of 
the brightest regions in the HERITAGE SDP strip. Figure [2] 
is a three-color composite image of N113 combining SPIRE 
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Fig. 2. Three-color composite image of the Hii region N113 
combining SPIRE 250 yum (red), PACS 160 fim (green), and 
MIPS 24 yum images. Our reliable Herschel sources are marked 
with red boxes. 5/7/^z^r-identified YSO candidates, YSOs "con- 
firmed" via spectroscopy or associated with OH and H2O masers 
("confirmed " YSOs) are also indicated. Maser p ositions are from 
iGreen et al.1 (.2008.) and EUingsen et al.l (|201qI) . 
SEDs of YSO-1 and YSO-2 in SectS 
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250 yum, PACS 160 yum, and MIPS 24 yum images. We detect 
9 reliable Herschel sources in N113, including 4 not previ- 
ously identified as YSO candidates. N 113 presents a clear ex- 
ampl e of star formation triggered by the winds from massive 
stars (lOliveira et al .120061) . Current star formation activity is con- 
centrated in the central lane of remnant molecular gas and dust, 
compressed by a complex structure of ionized gas bubbles, cre- 
ated by massive stars in several young clusters (< 10 Myr). 
Ongoing star formation in the lane is pinpointed in its earliest 
stages by maser e mission and by continuum emission as mas - 
sive YSOs evolve (iBrooks & Whiteoakll997l:IWong et al.l2006h . 
Two bright Herschel sources associated with masers and a third, 
fainter source between them are all confirmed YSOs. We model 
the full SED (1.25-500 yum) of the westernmost of these (YSO- 
1), to show HerscheVs impact on parameters of a known YSO. 
We also model a newly identified Herschel source YSO-2 that 
lies in a filament ~4' to the north (Fig.O online only). 

We present detailed analysis of two other sources, YSO-3 
and YSO-4. YSO-3, which appears to lie in a dense knot at 
the rim of a cavity in the molecular cloud of N 105, is asso- 
ciated with 6.67 GHz and 12.2 GHz methanol maser emission 
(Fig.m online only; Sinclair et al. 1992; Ellingsen et al. 2010). 
Methanol masers are closely associated with the earliest stages 
of massive star formation and are powerful probes of young mas- 
sive stars and protostars. We expect the SED shape of YSO-3 to 
be representative of SEDs of the youngest sources. YSO-4 is a 
new Herschel YSO candidate. It lies approximately 2' from the 
nearest Spitzer YSO candidates on the edge of an evacuated su- 
pershell in the LMC bar. The source is very faint shortward of 
IRAC 5.8 yum and then brightens in a "pillar of creation" mor- 
phology, indicative of a highly embedded YSO (Fig. [3 online 
only). 



5. Modeled physical properties 

Far-IR measurements are critical for deriving the total mass of 
circum(proto) stellar dust associated with a YSO, and inferring 
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Fig. 6. SEDs of four typical YSOs. Cyan diamonds are the flux 
values from Table [T] Cyan triangles indicate flux upper limits. 
The range of well-fitting models from the R06 YSO grid are 
shown in grey with a black line for the single best-fitting model. 
The best-fitting R06 model without Herschel {A <50 yum) is a 
dashed black line. A single-temperature greybody is shown fit 
long- wavelengths (green line), and a manually tailored radia- 
tive transfer model is shown as a red line. Model parameters are 
listed in Tabled 



its evolutionary state. Figure [6] shows typical sources that illus- 
trate our key findings (photometric data are listed in Table [T] 
online only). We fit each source with the R06 YSO model grid, 
first using only pre-if^r^c/z^/ photometry (/I < 50yum), and again 
using all data {A < 500 yum). Physical parameters derived from 
these fits are listed in Table [21 The addition of Herschel data 
provides a much tighter constraint on SED fits, and the range of 
well-fit models and uncertainties of derived parameters (as listed 
in the table) are thus significantly reduced. 

Nearly all sources detected by Herschel are Stage O/I, i.e., 
very young. A protostar's evolutionary state is often quantified 
by the ratio of circumstellar to stellar mass. In our single-source 
YSO models, the accreting flattened envelope has an analytic 
mass distribution (derived by Tereb ey et al. 1984) , and infall 
rates are proportional to envelope mass within an enclosed ra- 
dius. One caveat is that we do not know if all mass enclosed by 
the Herschel beam will fall on to the central source, so the de- 
rived infall rates may be upper limits. We find that the accretion 
rate relative to the cent ral mass, Menv/Mj, is a good measure of 
evolutionary state (e.g. lWhitney et al.ll2003al: llndebetouw et al.l 
2007; R06; R07). Without Herschel data, Menv/M^ of our 
sources is (2,5,l,4)xl0-5 yr'^ for (YSO-l,YSO-2,YSO-3,YSO- 
4). When Herschel photometry is included, Menv/^i^ increases 
to (2,3,2,4)xl0"^ yr"^ The higher infall rates are i n agreement 
with previous outflow studies of massive protostars (IChurchwelll 
Q999). The implied ages of <10'^ yr further strengthen the asser- 
tion that Herschel probes the very youngest sources. 

Observationally, the ratio of far-IR/submm luminosity to 
bolometric luminosity is often used as a proxy for the circum- 
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Table 2. Modeled YSO parameters, fitting each source using 
the R06 YSO grid and only pve-Herschel data (A < 50 yum), the 
R06 grid and all data (/I < 500 yum), and optically thin single- 
temperature dust emission using data A > 100 yum (greybody). 



Source Model 



A 



\og(Lboi) 

(Lo) 



log(M,^ 

(Mq yr- 







log(Mj,-,^) 

(Mo) 



M, 
(Mo) 



YSO-1 


R06 


<50 


5.1+0.2 


-3.3+0.8 


-1.0+1.1 


27+5 




R06 


<500 


5.5+0.1 


-2.1+0.0 


-0.1+0.4 


45+5 


YSO-2 


R06 


<50 


4.5+0.5 


-3.2+0.9 


-0.9+0.7 


13+2 




R06 


<500 


4.5+0.3 


-2.4+0.2 


-0.7+0.7 


13+2 


YSO-3 


R06 


<50 


4.2+0.2 


-3.8+1.8 


-1.2+1.1 


12+1 




R06 


<500 


3.9+0.2 


-2.5+0.2 


-0.8+0.8 


13+1 


YSO-4 


R06 


<50 


3.3+0.5 


-3.6+1.1 


-1.4+0.9 


7+1 



R06 


<500 


3.3+0.2 


-2.5+0.1 


-0.8+0.7 9+1 


Source Model 


A 


\0g(Lfar-IR) 

(Lo) 


^dust 

(K) 


log(M^^,) 

(Mo) 


YSO-1 greybody 


>100 


5.0 


20 


4.7 


YSO-2 greybody 


>100 


3.9 


19 


3.8 


YSO-3 greybody 


>100 


3.6 


24 


2.9 


YSO-4 greybody 


>100 


3.0 


16 


3.3 



stellar to total mass ratio. lAndre et aP (1 19931) define Class as 
having L^>350//m/L^o/ >5xl0"^. We calculate a similarly simple 
evolutionary measure, fitting an optically thin greybody to the 
data A >100yum, and find that it agrees reasonably with the YSO 
model fit results (Model "greybody" in Table [2] and green line 
in Fig.O. Resultant dust mass and far-IR/submm luminosity de- 
pend on the dust opacity power-law index. We use/5 = 2; a diff'er- 
ent choice does not change the conclusion that these sources are 
highly embedded. We find Lfar-m/l^boi ~ 0. 2 - 0.5 for our four 
examples, very high and consistent with the lAndre et al.l (1 19931) 
Class definition. We also find T^ust ~ 16-24 K, consistent with 
the far-IR colors (Fig. [T]) and expectations for the least evolved 
YSOs. 

Even the best- fitting R06 models (dashed line in Fig. [6]) often 
fall well below the data at longer wavelengths (also apparent in 
Fig. [B. The R06 models were designed for analysis of Spitzer 
data at shorter wavelengths; models have outer envelope radius 
^out ^ 0-5 pc, T> 30 K for luminous sources. The envelope of 
cooler dust and gas farther from the source emits little at /I < 100 
yum but can contribute significantly in the Herschel bands. 

We created new dust radiative transfer models (red line on 
Fig.O and find that implied accretion rate and other critical evo- 
lutionary parameters do not change significantly from the R06 
best fits. The new models have larger envelope radii, heating by 
external radiation, and R\H emission. We used the R06 best-fit 
models as starting points. Rout was varied between the R06 maxi- 
mum of 0.5 pc and 5 pc, with best-fitting values of (3.4, 1 , 1 , 1) pc. 
The importance of external heating wa s assessed by co mparing 
models with none to models with a lMathis et al.l (Il983h diff'use 
interstellar radiation field extincted by Av= 1 . In all cases the lat- 
ter improved the fits. Any increase in the outer envelope radius 
also increases dust absorption at short wavelengths, the eff'ect of 
which was balanced by varying the density in the outflow cavity 
between 0.1 and 1 times the R06 best-fit value. Reduced val- 
ues of (l,0.1,2,2)xl0"^^g cm"^ provided the best fits to the data. 
The envelope infall rates and source luminosities remained the 
same, showing that the R06 grid can still be used to interpret 
the Herschel data points. Adding R\H emission to our models 
results in better fits to bands eff'ected by aromatic emission (pre- 
dominantly IRAC 5.8 and 8.0 yum). Additional work must also be 



done to account for blending of multiple sources at Magellanic 
distances, as YSOs are frequently part of multiple systems and 
small unresolved clusters. 

Optimal analysis of Magellanic YSOs observed by Herschel 
and Spitzer will benefit from improvements to the R06 model 
grid. However, it is the inclusion of Herschel data, even with the 
older models, that is most significant: If these 4 embedded YSOs 
are typical, as we expect, the statistics of YSO physical parame- 
ters in the Magellanic System will be significantly clarified. 

6. Conclusions 

We show that Herschel far-IR photometry is very eff'ective in 
identifying YSOs in the LMC. Adding Herschel data to exist- 
ing Spitzer and near-IR observations results in significantly im- 
proved analysis of YSOs, as summarized in these key findings: 
o Nearly all sources detected by Herschel are Stage O/I, very 

young, with a high ratio of circumstellar to stellar mass. 
o Previously studied warm sources such as YSO-1 require 
more circumstellar dust to fit the Herschel data, implying a 
less evolved state than would be inferred from Spitzer data 
alone. 
o Herschel photometry significantly constrains our SED fits, 
decreasing the range of circumstellar dust masses and evolu- 
tionary states consistent with the measurements. 
o Many sources require even more cold circumstellar dust than 
is present in our original model grid, motivating improve- 
ments to our models. 
Our observations prove that Herschel can off'er us, for the first 
time, an inventory of the earliest stages of protostellar formation 
throughout an entire galaxy. 
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Fig. 3. Three-color composite images showing the environment of sources YSO-1 and YSO-2 in the N 1 13 massive star formation 
region based on the Spitzer and Herschel observations. The left-hand image combines IRAC and MIPS, showing emission from 
MIPS 24 yum, IRAC 8.0 yum, and IRAC 4.5 //m (red, green, and blue, respectively). The image at the right combines SPIRE 250 jim 
(red), PACS 160 yum (green), and MIPS 24 /^ m (blue) im ages. All images are scaled logarithmically. The linear distance scale 
is shown for the LMC distance of 50 kpc ( Schaefej 12008 ). Several categories of confirmed YSOs (see text for explanation and 
references) are marked as indicated in the legend. The OH and H2O masers marked by a single cross (x) symbol coincide in 
posi tion within the uncertainties. Sources YSO-1 and YSO-2 are discussed in Sec ts. |4] and [H Maser po sitions are from lGreen et al.l 
?2008) andEllingsenetal.'('2010V YSO-1 was identified a s a YSO candidate by Gruendl &Chu' (120091, 05 13 17.69-692225.0) based 
on the Spitzer data and confirmed spectroscopically by iSeale et al.l (i2009i) , although no ice features were detected. YSO-2 is a new 
Herschel YSO candidate. 
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Fig. 4. Three-color composite images showing the environment of YSO-3 in the LMC N 105 massive star formation region based 
on the Spitzer and Hersch el observations. Coloring and scaling are the same as in Fig. [3] YSO-3 was identified as a YSO candidate 
bv lGruendl&Ch u' ("20091 050958.52-685435.5) based on t he Spitzer data. This source is associated with 6.7 GHz and 12.2 GHz 
methanol masers (S inclair et al.lll992l:lEllingsen et al.ll2010h . 
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Fig. 5. Three-color composite images showing the environment of YSO-4 based on the Spitzer and Herschel observations. Coloring 
and scaling are the same as in Fig.O YSO-4 is located at the edge of an evacuated supershell in the LMC bar, between N 105 and 
N 1 13. It does not appear as a resolved source in 2MASS imaging and is faint in IRAC bands, excluding it from Spitzer YSO lists 
requiring shorter wavelengths. The SED of this cyidcnt Herschel source is quite steep from 24 jim to Herschel wavelengths (Fig.O, 
a trend we note in ~ 20% of the newly-identified Herschel YSO candidates 



Table 1. Spitzer and Herschel fluxes for YSOs discussed in Sect.|5] 



Source 


(hm s) 


Dec 
















Fluxes 


inmJy 










Fj 


Fh 


Fk 


^3.6//m 


FA.5nm 


Fs.Sjum 


Fs.Oium 


FMIPSlAjum 


FmiPSlQum 


FpACSimum 


FpACS\6Q^m 


FsPIRElSOjum FsPiE 


YSO-1 
YSO-2 
YSO-3 
YSO-4 


05:13:17.69 
05:13:17.99 
05:09:58.52 
05:13:01.55 


-69:22:25.0 
-69:20:04.9 
-68:54:35.5 
-69:10:34.8 


1.5(0.5) 

0.5(0.1) 

<0.2 

<0.1 


1.9(0.6) 

0.9(0.2) 

<0.5 

<0.1 


4.8(0.8) 

1.2(0.4) 

0.8(0.1) 

<0.3 


17(2) 

3.2(0.7) 

1.9(0.5) 

0.14(0.03) 


29(3) 

2.8(0.3) 

2.3(0.7) 

0.18(0.01) 


120(12) 
<24 
7(1) 
<1 


354(35) 

<68 

21(2) 

<3 


<3400 

580(96) 

250(88) 

8(1) 


26000(2600) 

2100(750) 

1500(680) 

310(73) 


53000(9600) 

3500(850) 

2100(300) 

400(51) 


66000(11000) 

5000(3500) 

1400(300) 

780(160) 


34000(3400) 1400C 

2800(460) 130C 

770(110) 430i 

630(63) 31C 



(a) Positions are determined from the highest-resolution images available (IRAC); few-arcsecond offsets present in this early Herschel data especially at PACSIOO are accounted for in the 

photometry. 

References. (D iGruendl &"Chul ( [2009h : (2) ISeale et al.l([2009h . 
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